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Abstract. AA Dor is an eclipsing, post common-envelope binary with an sdOB-type 
primary and an unseen low-mass secondary, believed to be a brown dwarf. Eleven years 
ago, a NLTE spectral analysis of the primary showed a discrepancy with the surface 
gravity that was derived by radial-velocity and light-curve analysis that could not be 
explained. 

Since then, emission Unes of the secondary were identified in optical spectra and 
its orbital-velocity amplitude was measured. Thus, the masses of both components are 
known, however, within relatively large error ranges. The secondary's mass was found 
to be around the stellar hydrogen-burning mass limit and, thus, it may be a brown dwarf 
or a late M-type dwarf. In addition, a precise determination of the primary's rotational 
velocity showed recently, that it rotates at about 65 % of bound rotation - much slower 
than previously assumed. 

A new spectral analysis by means of metal-line blanketed, state-of-the-art, non- 
LTE model atmospheres solves the so-called gravity problem in AADor - our result 
for the surface gravity is, within the error limits, in agreement with the value from 
light-curve analysis. 

We present details of our recent investigations on AA Dor 



1. Introduction 



AADor is a close, eclipsing binary ([Kilkenny et al.|[l978 ) and, thus, it is possible to 



determine its geometrical parameters (inclination /, radii and rggc) precisely by high- 
speed photometry. This has been done now over a time-span of about 50 000 eclipses 
and orbital period P = 0.261539 7363 (4)d (Kilkenny 2011) and / = 89?21 ± 0?30 
(iHilditch et al.ll2003h were determined. 



An early model of AADor ( Paczynskilll980h based on a light-curve analysis as- 



sumed an sdOB-type primary with Mpn = 0.36 Mq and a low-mass secondary with 
Msec - 0.054 Mq. Paczyriski estimated that the common envelope was ejected 5 • 10^ 
years ago and the primary will become a degenerate, hot white dwarf in another 5-10^ 
years. Within 5 • lO'" years then, gravitational radiation will reduce the orbital period 
to about 38 minutes and the degenerate secondary will overflow its Roche lobe, making 
AA Dor a cataclysmic variable. However, since AA Dor is a PCEB and is now a pre- 
cataclysmic variable, it is a key object for understanding such an evolution. Moreover, 
the low-mass companion, just at the limit of the hydrogen-burning mass, and its angu- 
lar momenti im are important for the understanding of the common-envelope ejection 



mechanism (iLivio & Sokerlll984l) . 
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Paczyiiskl (ll98Clh already pointed out the importance to find signatures of the ir- 
radiated secondary to m easure its radial-velocit y curve in order to refine the model of 
AA Dor. Quite recently, IVuckovic et all (l2008l) found spectral lines of the secondary 
in the UVESQ spectra (Progid 66.D-180) and measured its orbital-velocity amplitude 
(^sec > 230 km/sec), both components' masses are known now Mpn = 0.45 and 
Msec - 0.076 Mq, albeit with large error bars. 



Based on their spectral analysis of the sdOB primary, iKlepp & RauchI (1201 lb 

= 0.0725 - 0.0863 Mq within small er- 



found Mpri = 0.4714 



0.0050 Mq and M,e., 



ror limits. The scenario dXrimble & Aschwandenll200l1) 



"that a planet belonging to AA Dor tried to swallow its star during common- 
envelope binary evolution, rather than the converse (.Rauch..200ft) " 



(Msec - 0.066 Mq ~ 70 Ma.) is n ot valid anymore. 

The results of iKlepp & Rau ch (2011) for the primary, Teff 
log g - 5 .46 ± 0.05, are now in good agreement with the photometric model of lHilditch et al, 



: 42000 ± 1000 K and 



Detailed information about AA Dor and previ ous an alyses are given in iRauchI 

HOOO, 2004), Fleig et al. (2008), and Klepp & R auchI (l201ll). 

We will not repeat the description of the iKlepp & RauchI (1201 ih analysis here but 
we will emphasize problems that we encountered during that analysis. 



2. Atmosphere Modeling and Spectral Analysis of sd(0)B-type stars 



Miiller et"aLl ( 2010l ) presented a spectral analysis of the primary o f AAD or that resulted 



in a too-low T^fi. This initiated the re-analysis of iKlepp & RauchI (l201lh . They used the 
Tubingen Model- Atmosphere Package (rMA/|3) to calculate plane-parallel NLTE mod- 
els in hydrostatic and radiative equilibrium. Fig. [T] shows that the primary of AADor 
is located in the domain of static NLTE models. The limit between NLTE and LTE 
shown in Fig.[T]is not stringent. There are, however, deviations from LTE in any star, at 
least at high energies and high resolution, (cf. Jeffery, Pereira, Naslim, & Behara these 
proceedings.) 

In addition to the necessity to use NLTE models for an appropriate analysis of the 
primary of AADor, fully metal-line blanket ed model atmospheres should be used t o 
avoid the so-called Balmer-line problem (cf. iNapiwotzki & Rauchlll994l : lRauchll2000h . 



H4-He and the intermediate-mass elements C-i-N-i-O-i-Mg-i-Si-i-P-i-S were considered us- 
ing classical model atoms, while for Ca-i-S c-HTi-i-V-l-Cr-l-Mn-i-Fe-i-C o-i-Ni, the so-called 
iron-group elements, a statistical approach dRauch & Deetienl l2003 ) is employed to cre- 



ate model atoms. The most recent atomic data were used. 530 levels are 

treated in NLTE with 7 71 individual lines (from H - S) and 19 957 605 fines of Ca - Ni 
from Kurucz' line lists ( Kuruczir2009h combined to 636 superlines. For all models, we 



use the same model atoms and the same frequency grid. 

Two model-atmosphere grids were calculated, firstly a coarse grid (638 models) 
within Teff - 35 000 - 49 000 K (A T^s = 500 K) and log g - 5. 15 - 6.20 (A log g = 0.05). 



'Ultraviolet and Visual Echelle Spectrograph at ESO's VLT KUEYEN 

-http://astro.uni-tuebingen.de/~TMAP 

'http://astro.uni-tuebingen.de/~TMAD 
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Figure 1. Position of the primary of AA Dor in the log^ - Teff plane. Its symbol's 
size represents the errors of lKlepp & RauchI (1201 li) . The /u - I and jj - 4 lines indi- 
cate the Eddington limits for pure hydrogen and helium atmospheres, respectively. 



Secondly, a finer sub-grid (527 models) with Teff = 39 500 - 43 500 K (A Tss = 250 K) 
and Iog^ = 5.30 - 5.60 (Alog^ = 0.01). All calculations were performed on computa- 
tional resources of the bwGRiE0- 

The computational time for single models to converge is relatively long, and it took 
more than a diploma student's time to complete the grids. During their calculation, we 
performed quality control to check the grids' status. Fig. [2] shows a comparison of 
actual and desired T^s for two selected phases of model-atmosphere calculation. The 
left panel shows the situation just after a grid extension towards higher T^s and log g. 
The right panel shows a later phase where the complete grid (but the highest log g) is 
more homogeneously converged. 

It turned out that the temperature structures of the models (Fig.[3ll suggest that the 
models are converged much earlier than according to the T^s criterion. However, even 
at the converged state, a close look to the temperature structures of the models reveals 
a surprise. Fig.|4]shows that there is obviously some grouping in the temperature struc- 
tures within -2.9 < log m < - 1.9. The reason is that the formation depths of absorption 
edges and lines are dependent on T^s and log g. This is demonstrated for log g in Fig.|4l 
E.g. the absorption edge of H i « = 9 forms at log m w -2.05 in the log g = 5.30 model 
and at log m « -2.27 in the log ^ = 5.60 model. (The individual edge or line that is re- 



"^bwGRlD (http://www.bw-grid.de), member of the German D-Grid initiative, funded by the Ministry for 
Education and Research (Bundesministerium fiir Bildung und Forschung) and the Ministry for Science, 
Research and Arts Baden- Wiirttemberg (Ministerium fiir Wissenschaft, Forschung und Kunst Baden- 
Wiirttemberg) 
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Figure 2. log(rj™'*''' actual^ j-jnodei target^ ^.^^ points in time during the model-grid 
calculation. 



sponsible for the grouping is not investigated.) Our TMAP model atmospheres consider 
by default 90 depth points. In a test calculation, we used 360 depth points (Fig.© and 
could clearly show that the depth-point discretization is the reason of the temperature 
grouping. 

A higher number of depth points increases the calculation time linearly but has 
no significant impact on the analysis. However, one has to be aware that there may 
be artificial temperature jumps of a few hundred K in the line-forming regions of the 
atmosphere due to numerical limitations in the model-atmosphere calculations. This 
is especially important in case that fits are presented, with the extraordinary small 
statistical errors. 



3. The Balmer-Line Problem - Still present? 



The final Balmer-line fit (Fig. [5]) of lKlepp & RauchI ( 2011 ), got a comment by the A&A 
editor Ralf Napiwotzki: 

"It is slightly depressing that after all this improvements and on this high 
level of sophistication, the Balmer lines are still not well fitted." 



Indeed, the line cores as well as inner line wings of HyS - 5 are not well reproduced. 
Their line centers form at -2.4 < logm< - 1.8 (Fig. O where the temperature jumps 
(Fig. |4ll occur. A higher number of depths points, however, did not improve the agree- 
ment. 

The reason is more likely due to the Balmer-line problem (iNapiwotzki & Rauch 



19941) . i.e. that opacities are still missing in our model-atmosphere calculations. These 




Figure 3. Temperature structure of our retf = 42 000K models for all 
log g = 5.30 -5.60 (Alogg = 0.01). (logg = 5.30 is blue and logg = 5.60 is red in 
the online version.) The depth points (formation depths) where the line centers of 
selected hydrogen and helium lines become optically thin (r < 1) in the outer atmo- 
sphere are marked. 



are on the one hand of chemical species that are as yet unconsidered and on the other 
hand of species whose abundance is determined too-low. 

A new grid of models is already calculating, it will be used to further investigate 
this issue. 



4. Conclusions 



The sdOB primary of AA Dor is well studied and most of its parameters are determined 
within small error ranges. Additional opacities may improve the model atmospheres 
and e.g. the agreement of synthetic Balmer-line profiles with the observation. We 
summarize parameters of A A Dor in Tableland show a sketch of the binary in Fig.[6l 
The masses of both components of AADor are well known, Mpn = 0.4714 ± 
0.0050 Mo and Msec = 0.0725-0.0863 Mq (Klepp & Rauc hl201 ll). Since the hydr ogen- 
burning mass limit is ~ 0.075 Mq (iChabrier & Baraffelll997l : IChabrier et al.ll200d) . it is 
still unclear whether the secondary is a brown dwarf or a late M-type dwarf. In order 
to make progress and to investigate on the nature of the secondary of AA Dor, phase- 
dependent spectra of AA Dor in the infrared are highly desirable in order to identify the 
secondary's spectral contribution and to follow its spectral evolution during its orbit. 
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Figure 4. Top panel: Detail of Fig. [3] The formation depths of the hydrogen ab- 
sorption edges are marked at top (for log ^ = 5.30, blue) and bottom (log ^ = 5.60, 
red). Bottom panel: Same as top panel but with a four times refined depth-point 
discretization in the model. 



5. An Answer to an Important Question 

At the very end of the question time following my talk, Darragh O'Donoghue asked 

"Thomas, is it right, that with all your spectral analysis, you arrived at the 
same result that we got already 30 years ago from light-curve analysis?" 

Of cours e, he rneant t he SAAO group and log g of the primary - and thus, the log g 
problem (Rauch 2000), that is solved now. The simple answer to Darragh's question is 
yes. However, I wo uld like to poin t out that - thanks to the comments of Ron Hilditch, 
who was referee of iRauchI ( 20001) - we took the challenge to solve the logg problem 



and finally achieved agreement. This is a good example that, if carried out with care, 
theory works. The independent method to determine logg from light-curve analysis 
provided a crucial constraint to test our model atmospheres. Within this course, we 
improved o ur models, the atomic data, as well as our observational, spectral data base 
(lRauch& Werner. ,2003) . The errors in Teff, logg, and abundances (iFleig et al.l 120081) 
for the primary of AA Dor were reduced. 

In total, our picture of AA Dor improved, making this one of the best-analyzed 
PCEB. We are still not tired to improve the primary's model, but it appears time to 
focus now on the outstanding question of the nature of the low-mass companion. Phase- 
dependent, high-S/N infrared spectroscopy, that covers the complete orbital period of 
A A Dor, will provide answers. 
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Figure 5. Comparison of theoretical line profiles of H and He lines calculated 
from the fi nal Teff =42000 K and log ^ = 5.46 model with the UVES observation (for 
details, see iKlepp & RaucEIMlh . 
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Figure 6. Dimension of AA Dor. 
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Table 1. P arameters of AA Dor compiled f romlRauchI (l2000t) . lFleig et alj (l2008h . 
iMuller et all ^10), and iKlepp & RauchI (1201 ilk log ex are normalized to 
logZ/^xSx = 12.15. [X] denotes log(abundance/solar abundance) of species X. 
The solar abundances are given by Asplund et al. (2009). 



Primary 



Teff/kK 42.0 ±1 
logig/f) 5.46 -0.05 





mass 


number 






abundances 


fraction 


fraction 


log fix 


[X] 


H 


9.94E-01 


9.99E-01 


12.144 


0.130 


He 


2.69E-03 


6.80E-04 


8.977 


-1.967 


C 


1.78E-05 


1.50E-06 


6.320 


-2.124 


IN 


A 1 ^'^ 

'f . 1 jii-Uj 


a AAC AA 


D.DZl 


1 in. 

— 1 .ZZj 


o 


l.OlE-03 


6.39E-05 


7.950 


-0.754 


Mg 


4.08E-04 


1.70E-05 


7.375 


-0.240 


Si 


3.05E-04 


l.lOE-05 


7.186 


-0.339 


P 


5.20E-06 


1.70E-07 


5.375 


-0.050 


S 


3.24E-06 


1.02E-07 


5.155 


-1.980 


Ca 


5.99E-05 


1.52E-06 


6.325 


-0.030 


Sc 


3.69E-08 


8.33E-10 


3.065 


-0.100 


Ti 


2.78E-06 


5.89E-08 


4.915 


-0.050 


V 


3.73E-07 


7.42E-09 


4.015 


0.070 


Cr 


1.66E-05 


3.24E-07 


5.655 


0.001 


Mn 


9.88E-06 


1.82E-07 


5.405 


-0.040 


Fe 


1.15E-03 


2.09E-05 


7.465 


-0.050 


Co 


3.59E-06 


6.17E-08 


4.935 


-0.069 


Ni 


3.48E-04 


6.01E-06 


6.923 


0.689 


Mpri/Mo 


0.4714 


±0.0050 






■'^pri / R© 


0.21 


+0.028 
-0.024 






Opri / R© 


0.208 


+0.030 
-0.018 






^pri / Lq 


120 


+ 15 

-20 






^rot/ sec 


30 


±1 







Secondary 



M,ec/Mo 0.0788 ^^^"^^^ 

Rsec/^Q 0.11 -aoi3 

«sec/R© 1.244 ^^aiio 

d/pc 352 +20 

«H/cm2 2-1020 ±1-1020 

^B-v 0.01 ±0.01 
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David Kilkenny, with small tribal face painting, at Spier wine estate, July 27, 2011. 



sdOB5 marked the 65 birthday of David Kilkenny, who discovered that AA Dor is a 
short-period, eclipsing binary system. Since December 1974, he observed this system 
once in a while, covering now a period of more than 46 119 eclipses, and could show 
that the period increase or decrease is less that 10"'^ d/orbit. 

All the best wishes to you, Dave, and have the energy for the next tens of thousands 
of eclipses of AADor! 



